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ABSTRACT 

We report the discovery of a cluster-scale lensed quasar, SDSS J1029+2623, selected from the Sloan 
Digital Sky Survey. The lens system exhibits two lensed images of a quasar at z s — 2.197. The image 
separation of 22'.'5 makes it the largest separation lensed quasar discovered to date. The similarity of 
the optical spectra and the radio loudnesses of the two components support the lensing hypothesis. 
Images of the field show a cluster of galaxies at zi ~ 0.55 that is responsible for the large image 
separation. The lensed images and the cluster light center are not collinear, which implies that the 
lensing cluster has a complex structure. 

Subject headings: gravitational lensing — quasars: individual (SDSS 102913.94+262317.9) — galaxies: 
clusters: general 
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1. INTRODUCTION 

The discovery of SDSS J1004+4112 with an image 
separation of 14'.'6, the first example of a quasar mul- 
tiply imaged by a massive cluster of galaxies, opened a 
new window for understanding our universe (Inada et 
al. 2003; Oguri et al. 2004; Sharon et al. 2005). Al- 
though there are many examples of galaxies (arcs) lensed 
by clusters, large-separation lensed quasars have several 
advantages over arcs as a cosmological probe. First, the 
simpler (point-like) structure of quasars and their well- 
understood redshift distribution should make the large- 
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separation lensed quasars much cleaner probes of cos- 
mology and structure formation models (e.g., Oguri & 
Keeton 2004; Hennawi et al. 2006), while the statistics 
of arcs remain contentious (Bartelmann et al. 2003). Sec- 
ond, the time-variability of quasars allows the measure- 
ment of time delays among the multiple lensed images, 
thereby breaking the mass-sheet degeneracy of lens mod- 
els given a priori knowledge of the Hubble constant (e.g., 
Kochanek 2002). This was explored in detail for SDSS 
J1004+4112 (e.g., Oguri et al. 2004; Williams & Saha 
2004; Kawano & Oguri 2006; Fohlmeister et al. 2006). 
Our current problem is that the small number of known 
systems limits their use in statistical analyses, however, 
current and future large lens surveys will discover many 
large-separation lensed quasars (Wambsganss 2003). 

In this Letter, we report the discovery of the sec- 
ond large-separation lensed quasar, SDSS J1029+2623, a 
quasar at z s = 2.197 doubly imaged by a massive galaxy 
cluster at z;~0.55. It was discovered in the course of the 
Sloan Digital Sky Survey Quasar Lens Search (SQLS; 
Oguri et al. 2006), which is a survey of strongly lensed 
quasars in the Sloan Digital Sky Survey (SDSS; York et 
al. 2000). The image separation of 22'.'5 makes it the 
largest separation lens among the ~ 100 lensed quasars 
known so far (Kochanek 2006). 

2. CANDIDATE SELECTION AND SUBARU 
OBSERVATIONS 

From the SDSS-II Sloan Legacy Survey data, we dis- 
covered the lens candidate SDSS J1029+2623 using the 
algorithm described in Oguri et al. (2006). It consists of 
a quasar whose redshift was spectroscopically measured 
by the SDSS to be z = 2.198 (hereafter referred to as 
A) and a stellar object (hereafter B) with almost the 
same color as component A (see Table 1). Components 
A and B are separated by 22'.' 5. There are no other point 
sources in the field with similar colors that would be can- 
didates for additional quasar images. In addition to the 
similar colors, the gri composite SDSS image shown in 
the left panel of Figure 1 exhibits a concentration of red 
galaxies centered on two bright galaxies (named Gl and 
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Fig. 1. — Left: The gri composite SDSS image of SDSS J1029+2623 (1'.'3 seeing). The quasar images (blue stellar objects) are indicated 
by A and B. Gl and G2 (red extended objects) are likely to be member galaxies of a lensing cluster at z~0.55. The inset shows an expanded 
view of component B: An object ~ 2" Southeast of component B has quite different color from those of the quasar components. Middle: 
The FOCAS i-band image (l'.'S seeing). Right: Locations of galaxies (i < 22.5) identified in the FOCAS images are plotted with filled 
circles. Larger circles mean brighter galaxies. Galaxies that survive the color cut of R — i = 0.7±0.3 are shown by darker circles. Filled 
triangles are the locations of the two quasar images. In all the panels, North is up and East is left. 

G2), which implies the existence of a high-redshift clus- 
ter of galaxies. Indeed, the photometric redshifts of Gl 
and G2 are both ~0.55 from their SDSS colors (Csabai 
et al. 2003). The properties of components A, B, Gl, 
and G2 are summarized in Table 1. We note that the 
details of the SDSS are described in a series of technical 
papers: Gunn et al. (2006) for the dedicated wide-field 
(3° field of view) 2.5-m telescope; Fukugita et al. (1996), 
Gunn et al. (1998), Lupton et al. (1999), Hogg et al. 
(2001), Lupton et al. (2001), Smith et al. (2002), Pier et 
al. (2003), Ivezic et al. (2004), and Tucker et al. (2006) 
for the photometric survey; Richards et al. (2002) for 
the spectroscopic target selection algorithm of quasars; 
and Blanton et al. (2003) for the tiling algorithm of the 
spectroscopic survey. Most of the SDSS data are already 
publicly available (Stoughton et al. 2002; Abazajian et 
al. 2003, 2004, 2005; Adelman-McCarthy et al. 2006). 

We obtained 600 sec spectra of components A and B 
with the Faint Object Camera And Spectrograph (FO- 
CAS; Kashikawa et al. 2002) at the Subaru 8.2-meter 
telescope on 2006 June 28. The observation was con- 
ducted in the 3x1 on-chip binning mode, using the 300B 
grism, the SY47 filter and a O'.'8-width slit aligned along 
components A and B. The spectral resolution, wave- 
length coverage, and spatial scale of the CCD detector 
were R ~ 500, from 4700 A to 9400 A, and 0'.'311 pixel" 1 , 
respectively. Although the seeing was poor (~1'.'8), the 
22'.'5 separation makes it easy to extract the spectra of 
the two components using standard IRAF 18 tasks. The 
spectra are shown in the upper panel of Figure 2. It is 
clear that both components are quasars at the same red- 
shift, z = 2.197 (see Table 1). Moreover, the spectral 
shapes are similar; both components have similar broad 
absorption line features at the same wavelength in the 
C IV emission lines, similar profiles for the Fe III emis- 
sion lines, and similar characteristic shapes of the red 
wings of the Mg II emission lines. In addition, the spec- 
tral flux ratio plotted in the lower panel of Figure 2 is 
constant (~1.2) over the full range of the observed wave- 
lengths. The similarity of the spectra supports the idea 
that the two components are lensed images of a single 
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Fig. 2. — Spectra of components A (black solid line, shifted 
upwards by a constant of 10 flux density units) and B (gray solid 
line) taken with the FOCAS on the Subaru telescope. The signal 
to noise ratio (in the 6000 A - 8000 A continuum) is about 25.0 per 
pixel. The dotted lines indicate the positions of the quasar C IV, 
C III], Fe III, and Mg II emission lines redshifted to z = 2.197. 
The feature at 7600 A is atmospheric. The inset shows an expanded 
view of the C IV emission lines. The vertical dot-dashed line in the 
inset shows the central position of the broad absorption lines. In 
the bottom panel, the black solid line shows the spectral flux ratio 
between A and B, and the dark gray solid line shows the ltr error of 
the spectral flux ratio, which is derived from the noise per pixel in 
each spectrum. The ratio is confirmed to be constant (~1.2) within 
the 2.3(7 error, over the entire range of observed wavelengths. 



18 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



quasar. 

We also obtained 120 sec R- and 120 sec i-band im- 
ages (~l'/8 seeing) with FOCAS. On-chip 2x2 binning 
yielded images with a pixel scale of 0'.'207 pixel -1 . The 
i-band image is shown in the middle panel of Figure 1. 
From the FOCAS images, we found that galaxy Gl is 
likely to be a superposition of two early-type galaxies; 
however, we treated them as a single object since the 
poor seeing prevents us from separating them. We used 
SExtractor (Bertin & Arnouts 1996) for the photome- 
try (MAG_AUTO) calibrated by the SDSS magnitudes 
of the nearby stars 19 . All galaxies (classified by the com- 

19 We used the transformation formulae described in Table 7 of 
Smith et al. (2002) to calculate R magnitudes of the nearby stars 
from the SDSS gr magnitudes. 
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Fig. 3. — Constraints on the center of the potential of the lens- 
ing cluster from mass modeling. Solid lines indicate 1, 2, and 3<r 
likelihood regions for the center based on an SIE model for the 
cluster potential (see text for details). The quasar images (A and 
B) and bright member galaxies (Gl and G2) are marked by filled 
triangles and crosses, respectively. If we model the system using a 
SIE+external shear with the potential center fixed to galaxy Gl, 
the best-fitting model has an Einstein radius of 15"0 (correspond- 
ing to the velocity dispersion of ~900 km sec -1 ), an ellipticity 
of 0.25 with its direction of 88° (East of North), and a shear of 

0. 09 with its direction of 28°, although it does not reproduce the 
observables exactly. 

bination of the SExtractor CLASS-STAR parameter and 
visual inspection) with i < 22.5 are plotted as filled cir- 
cles in the right panel of Figure 1. Among them, we 
selected galaxies with colors of R — i = 0.7±0.3 as possi- 
ble cluster members, because Gl and G2 have R—i colors 
of 0.6 and 0.8, respectively. The R — i color of 0.7 is con- 
sistent with the red-sequence color of a galaxy cluster at 
z~0.55 (e.g., Goto et al. 2003). The selected galaxies are 
plotted as dark filled circles in the right panel of Figure 

1. As expected, clustering of galaxies near Gl and G2 
is much more pronounced after the color cut, suggesting 
the existence of a cluster of galaxies at z ~ 0.55. 

3. DISCUSSION AND CONCLUSION 

An interesting thing about this lens system is that the 
quasar images and the light center of the lensing cluster 
appear not to be collinear (see Figure 1). Brightest clus- 
ter galaxies are sometimes significantly offset from the 
center of a cluster mass (e.g., Lin & Mohr 2004), so the 
lack of collinearity does not necessarily argue against the 
lensing hypothesis. Nevertheless, it is of great interest 
to understand the impact of this offset for lens models. 
To explore this, we use standard mass modeling tech- 
niques as implemented in the public lensmodel software 
(Keeton 2001). Since the small number of observational 
constraints limits detailed investigations of the mass dis- 
tribution, we adopt the simplest model used for model- 
ing studies, a singular isothermal ellipsoid (SIE). We fit 
the quasar positions and fluxes assuming a position er- 
ror of 075 and a flux error of 20%; we adopt errors larger 
than the actual measurement errors in order to account 
for perturbations from massive galaxies in the cluster. 
In addition we include a mild prior on the ellipticity, 
e = 0.3 ± 0.2, to exclude models with a highly elongated 
cluster. With this model, we derive the allowed range 
for the center position of the lensing cluster potential, 
as shown in Figure 3. The offset of at least ~ 30 kpc 
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Fig. 4.— The 6 cm VLA radio map of SDSS J1029+2623 field. 
North is up and East is left. The observation was made during the 
transition between the B and CnB arrays, leading to the strongly 
elongated beam (shown in the lower left inset). The red small 
dots represent the optical position of components A, B, Gl, and 
G2. The small blue dot Northwest of G2 marks the ROSAT X-ray 
source position while the surrounding blue dashed circle indicates 
the uncertainties in the position. 

between the preferred model positions and galaxies Gl 
and G2 suggests either that the actual mass distribution 
of the lensing cluster is more complex than this simple 
model or that there is a genuine offset that might be 
confirmed in X-ray observations. 

While the similarity of the optical spectra and the exis- 
tence of the lensing galaxy cluster make the lensing inter- 
pretation solid, we can further test the lensing hypothesis 
using archival data at other wavelengths. In the FIRST 
radio survey (Becker et al. 1995), there are two radio 
sources in the field: Galaxy G2, with a peak flux of 4.20 
mjy/beam, is the strongest source, component B, with a 
peak flux of 1.01 mjy/beam, is detected just at the flux 
limit, and component A is not detected. Although the 
difference could be explained by quasar variability com- 
bined with the time delay between the lensed images, 
it might argue a binary interpretation for the quasars. 
In order to clarify the interpretation of the FIRST ra- 
dio data, we obtained a 1 hour 6 cm Very Large Array 
(VLA) radio map of the system on 23 September 2006. 
Figure 4 shows the primary beam corrected VLA 6 cm 
map; both quasar components have more than 5<r de- 
tections (0.237±0.045 mjy/beam for A and 0.325±0.045 
mjy/beam for B). There are some offsets (but within 5<7, 
see Table 1) between the SDSS positions and the VLA 
positions of the quasar components. However both com- 
ponents are extended in the VLA data and better quality 
data will be necessary to clarify the offsets. The radio 
loudnesses of the images, defined relative to the SDSS 
i-band fluxes, are 1.8 and 2.5 for images A and B, re- 
spectively. The radio flux ratio (A/B) of 0.73 is mildly 
inconsistent with the optical flux ratio of 1.2, however, 
allowing for variability both between the epochs of the 
two observations and the multi-year time delay of the two 
images, finding two quasars with such similar radio loud- 
nesses further supports the lensing hypothesis given the 
enormous range of radio loudnesses observed for quasars 
(Ivezic et al. 2002). 

The system was also detected as 1RXS 
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J102912.0+262338 by the ROSAT All-Sky Survey 
(Voges et al. 1999) with a count rate of 0.040 ± 0.016 
counts s -1 . This corresponds to an 0.1-2.4 keV flux of 
approximately (4 ± 1) x 10~ 13 ergs cm~ 2 s -1 , and it is 
flagged as an extended source, albeit at low significance. 
The X-ray position could be consistent with the position 
of G2 given the lO'.'O uncertainties in the X-ray data 
(Figure 4). Assuming the X-ray signals arise from 
the lensing cluster (2~0.55), the X-ray luminosity is 
approximately 3 x 10 M h~ 2 ergs s -1 (for a cosmological 
model with Qm = 0.3 and = 0.7) and corresponds 
to a velocity dispersion of approximately a ~ 1000 km 
s _1 that is slightly larger than the velocity scale of 
a ~ 900 km s _1 needed to produce the 22'.'5 image 
separation. Thus the X-ray luminosity is a little too 
high, suggesting that the AGN in galaxy G2 is an X-ray 
source as well as a radio source, even though it shows no 
signs of AGN activity in its broad band optical colors. 

In summary, the best interpretation of the two quasar 
components in the SDSS J1029+2623 system is that a 
quasar at z s = 2.197 is doubly imaged by a cluster of 
galaxies at zi ~ 0.55. The evidence for strong lens- 
ing comes from the remarkable similarity of the spec- 
tral shapes and the similar radio loudnesses of the two 
quasars, and the existence of the lensing cluster with 
the presence of an extended X-ray source which is capa- 
ble of producing the observed image separation. How- 
ever, simple mass models require a center of the cluster 
potential offset from the positions of the bright galax- 
ies or the ROSAT X-ray emission. Such an offset with 
respect to the bright galaxies was also found for the 
SDSS1004+4112 system (Oguri et al. 2004), and the X- 
ray position may be dominated by the AGN in galaxy G2. 
Thus, further follow-up observations, such as deep and 
high-resolution X-ray imaging and optical spectroscopy 
of member galaxies are needed to identify the cluster po- 
tential center. 
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TABLE 1 
Properties of SDSS J1029+2623 



Component 


R.A.(SDSS) 


Decl.(SDSS) 


R.A.(VLA) 


Decl.(VLA) 


i 


u — g 


g-r 


r — i 


i — z 


Redshift a 


A 


10:29:13.94 


+26:23:17.9 


10:29:14.05 b 


+26:23:17.4 b 


18.59+0.01 


0.72+0.04 


0.18+0.01 


0.19+0.01 


0.32+0.03 


2.1966+0.0003 


B 


10:29:14.24 


+26:23:40.1 


10:29:14.31 b 


+26:23:39.3 b 


18.61+0.01 


0.67+0.04 


0.14+0.01 


0.20+0.01 


-0.10+0.04 


2.1969+0.0003 


Gl 


10:29:13.35 


+26:23:32.8 






19.01+0.04 


0.71+0.90 


1.93+0.17 


1.05+0.06 


0.75+0.09 


0.56+0.03 


G2 


10:29:12.48 


+26:23:32.0 


10:29:12.49 


+26:23:31.9 


18.77+0.03 


1.63+1.72 


1.87+0.15 


1.06+0.05 


0.48+0.08 


0.53+0.06 



Note. — Data are from the SDSS except for the celestial coordinates in the VLA data and the spectroscopic redshifts. Magnitudes refer model magnitudes without Galactic extinction 
corrections, measured on 13 December 2004. All celestial coordinates are given in J2000. 



Spectroscopic redshifts (derived from Mg II emission lines in the FOCAS spectra) for A and B, and photometric redshifts for Gl and G2 (Csabai et al. 2003). 
b VLA position errors for A and B are both about 0'.'5 per coordinate. 



